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We describe traffic flows in one lane roadways using kinetic theory, with special emphasis on 
the role of quenched randomness in the velocity distributions. When passing is forbidden, growing 
clusters are formed behind slow cars and the cluster velocity distribution is governed by an exact 
Boltzmann equation which is linear and has an infinite memory. The distributions of the cluster 
size and the cluster velocity exhibit scaling behaviors, with exponents dominated solely by extremal 
characteristics of the intrinsic velocity distribution. When passing is allowed, the system approaches 
a steady state, whose nature is determined by a single dimensionless number, the ratio of the passing 
time to the collision time, the two time scales in the problem. The flow exhibits two regimes, a 
laminar flow regime, and a congested regime where large slow clusters dominate the flow. A phase 
transition separates these two regimes when only the next-to-leading car can pass. 



I. INTRODUCTION 

Traffic flows are strongly interacting many-body sys- 
tems. Therefore, theoretical techniques such as kinetic 
theory and hydrodynamics are useful in describing the 
rich phenomenology of traffic flows which includes shock 
waves, phase transitions, clustering, metastability, hys- 
teresis, etc. Traffic is typically modeled within macro- 
scopic descriptions such as hydrodynamics and kinetic 
theory &5|or microscopic approaches, e. g., c ellular au- 
tomata |3-12| and car- following models jKJ-|f7]]. The 
large body of recent work on the physics of traffic flows 
is surveyed in Refs. |I^-|20|l. 

In this review, we describe how quenched randomness 
in the car velocities leads to formation of clusters in one 
lane roadways. We assume ballistic motion with infinite 
memory, namely that each car has a preferred "intrinsic" 
velocity by which it drives in the absence of other cars 
ff] §6|. While the emerging behavior is quite similar to 
that found in stochastic particle hopping processes with 
quenched disorder p7|-p3[ , these ballistic motion models 
often admit deeper analytical treatment. Our starting 
point is an idealized no passing flow, where an exact an- 
alytical solution is possible, and an exact kinetic theory 
can be constructed |22|. We then treat more realistic 
generalizations where passing is allowed using approxi- 
mate kinetic theories [p3|-p6[. Our goal is to provide a 
concise summary where key features are emphasized and 
outstanding open issues are highlighted. 



II. NO PASSING ZONES 

Our basic traffic model mimics cluster formation (of- 
ten also called platoon formation) in one-lane roadways 
where passing is forbidden |^l|,^2| . In this model, each car 
moves ballistically at its initial velocity until it overtakes 
the preceding car or cluster. After this encounter, the in- 
cident car assumes the velocity of the cluster which it has 
just joined. Cars are taken to be size-less, and collisions 



to be instantaneous. We primarily consider spatially ho- 
mogeneous situations where the positions and the veloc- 
ities of the cars are initially uncorrelated. Specifically, 
cars are distributed randomly in space with a concen- 
tration Co, and their velocities are drawn from the initial 
velocity distribution Pq(v). Remarkably, analytic expres- 
sions can be obtained for the velocity distribution and the 
joint size-velocity distribution of clusters. Furthermore, 
it is also possible to describe analytically spatial inhomo- 
geneities and even input of cars into the system. 



A. The velocity distribution 



We first consider the cluster velocity distribution. In 
this description, only the lead car is relevant and trailing 
cars in a cluster can be ignored. Let P(v , t) be the distri- 
bution of clusters with velocity v at time t. Initially, all 
cars are lead cars and the cluster (or lead-car) velocity 
distribution equals P{v,t) — P (v)S(v,t), a product of 
the initial velocity distribution and the survival proba- 
bility S(v,t). The survival probability is the probability 
that a car with initial velocity v avoids "collisions" with 
slower cars up to time t, and hence, is still moving at the 
same velocity. Consider a car of initial velocity v. To 
ensure that it would not overtake slower cars of velocity 
v' < v, an interval of size (v — v')t ahead of it must not 
contain v'-cais initially. For the velocity distribution, 
Pq(v), and a Poissonian initial spatial distribution, the 
probability for such an event is exp [— t (v — v')Pq(v')]. 
For a car to survive to time t, this exclusion probability 
should be taken into account for every v' < v, and taking 
the product over all possible slower cars yields the sur- 
vival probability, S(v,t) = exp[— t j Q dv'(v — v')P (v')]. 
Consequently, the cluster velocity distribution is found 
for arbitrary initial conditions 



P(v,t) =P {v)exp 



t / dv'{v-v')P (v') 



(1) 
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The process is deterministic (the initial condition is the 
only source of randomness) and given the initial positions 
and velocities of the cars, the state of the system at any 
later time follows. This is reflected nicely in Eq. ([!]). 

The exact solution (0) satisfies the following Boltz- 
mann equation, 

^^- = -P(v,t)j V dv , (v-v')P(v',0). (2) 

Interestingly, this rate equation is linear in the velocity 
distribution P(v, t). The collision rate is proportional to 
the relative velocity, v — v' , and the initial velocity dis- 
tribution of the slower cars, a signature of the infinite 
memory in the system. This is a unique case where the 
hierarchy of evolution equations corresponding to the ve- 
locity distributions terminates at the first order. 

Whereas the steady state of one-lane traffic with no 
passing is trivial, viz. all cars will eventually join a cluster 
led by the slowest car in the system, the time dependent 
behavior is interesting. We concentrate on the long time 
behavior, which is largely independent of the initial dis- 
tribution of fast cars, as follows from Eq. (|l|). For discrete 
velocity distributions, the time dependence is exponen- 
tial, and we focus on continuous distributions. In this 
case, we find directly from Eq. (|l|) that both the cluster 
concentration, c(t) = J dvP(v,t), and the average veloc- 
ity, (v(t)) — cr 1 J dvv P(v,t), decay algebraically with 
time 



c (t) ~ r a 



H+l 
M + 2' 
1 

M + 2' 



(3) 



The scaling exponents a and (3 depend only on the small- 
v extremal statistics [[m] of the initial velocity distribu- 
tion, 



Po(v) 



0. 



(4) 



via the cutoff exponent fj, > — 1. The two scaling expo- 
nents are related by a + (3 = 1 as dictated by an elemen- 
tary mean free path argument: cvt ~ 1. 

Since the number of cars is conserved, the average clus- 
ter size is inversely proportional to the concentration, 
(m) ~ 1/c, and the size growth law is (m) ~ t a . In 
the limit fi — > oo, the size grows linearly with time. 
In contrast, when fx — * —1, the size remains roughly 
constant, since the velocity distribution becomes effec- 
tively unimodal and collisions become exceedingly rare. 
This qualitative dependence on the form of the initial ve- 
locity distribution is reminiscent of ballistic annihilation 
processes, where ballistically moving particles annihilate 
upon collision p5|-[3S(] . The above clustering process can 
be viewed as a ballistic aggregation process that possesses 
a single mass conservation law. The sensitive dependence 
on the initial conditions is in contrast with momentum 



conserving ballistic agglomeration processes (that mim- 
ics large scale formation of matter in the universe) where 
a universal scaling asymptotic behavior emerges [|9| [IlJ . 

The average velocity is the only relevant velocity scale 
in the problem and asymptotically the velocity distribu- 
tion follows the scaling form 

P(v,t) ~^- Q $(wt~' 3 ). (5) 

From Eqs. (|) and (||), the scaled distribution is &(z) = 
az^expi-bz^ 2 ), with b = o/[(/i + l)(/i + 2)]. There- 
fore, the small-w asymptotics of the initial velocity distri- 
butions governs not only the scaling exponents but also 
the entire shape (including the large velocity tail) of the 
scaling function $>(z). This scaling behavior indicates 
that at time i, most cars moving initially with veloci- 
ties larger than the typical velocity scale, (v(t)) ~ , 
have already joined clusters led by slower cars, while cars 
slower than this velocity scale are still driving with their 
initial velocity. 



B. The size- velocity distribution 

Given the nature of the model, a car is only affected 
by the initial configuration of cars ahead of it. This key 
feature enables solution of the cluster velocity distribu- 
tion and it allows treatment of a more detailed quantity, 
the joint size-velocity distribution. To obtain P m (v,t), 
the density of clusters of size m and velocity v, it is use- 
ful to introduce the cumulative distribution, Q m {v,t), 
the distribution of clusters of velocity v containing at 
least m cars. Knowledge of this cumulative distribution 
yields the joint size- velocity distribution via differencing, 
P m {v,t) = Qm(v,t) - Qm+l («>*)< 

Consider the first nontrivial quantity, Qa(v, £), the 
probability distribution of clusters of velocity v with at 
least two cars at time t. This quantity is equal to the 
product of the probability that lead car has survived up 
to time t, P(v,t), and the probability that the car trail- 
ing it actually experiences a collision prior to time t. Let 
x\ and v\ be the initial position and the initial veloc- 
ity of the trailing car, respectively. For such a collision 
to occur, the trailing car must be faster than the lead 
car, v\ > v, and the interval separating the two cars 
must be initially free of other cars. The probability for 
this composite event is the product of the probabilities 
of each individual event. Given a random (Poisson) spa- 
tial distribution, the probability an interval is empty is 
exponential in its length, and the collision probability is 



s (t>,t) =P(v,t) / duiPo(vi) 



xi<(vi—v)t 



dxi exp(— xi). 
(6) 



The fact that the trailing car cannot be slowed down by 
any other car before colliding with the lead car is crucial 
in obtaining this solution. This solution can be gener- 
alized to arbitrary cluster sizes. Following the two-car 
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cluster case, one simply integrates over all the initial po- 
sitions and velocities of the consecutive cars to eventually 
collide with the lead car in the cluster to give 



c(x,t) — t a C(X). Far from the origin, the scaled den- 



sity decays as C(X) ~ X implying c(x,t) 



for 



m ~ 1 /■OO p 

t)=p(v,t)n dviP (vi) / 



^1 H <(Vi —V)t 



x t a . Consequently, the total number of clusters in 

the originally empty region, N(t) = J °° dxc(x,t), grows 

logarithmically slow with time 
dxi exp(— Xi). 

Int. 



(7) 



The integration limits reflect the fact that all the colliding 
cars must be faster than the lead car, and the restriction 
on the integration limits ensures that cars are sufficiently 
close to the lead car so that collisions indeed occur. 

Given the cumulative car distribution, the joint size- 
velocity distribution can be formally obtained. Since 
(to) ~ t a and (v) ~ t _/3 , we anticipate the following 
scaling behavior: P m (v, t) ~ ^-^^(mf" , vt 13 ). This 
indeed holds and from Eq. (Q) one obtains the scaled 
joint distribution 



z) = cz^(x + zY l+1 exp \-b(x + zf 



(8) 



with c = a 2 / (fj, + 1). Again only two parameters, a and 
fi characterizing the small velocity characteristics of the 
initial conditions, are needed to fully describe the asymp- 
totic state of the system. The joint distribution (||) pro- 
vides a comprehensive description of the traffic clustering 
process. It may be considered as the counterpart of the 
well-known result for diffusion-controlled aggregation in 
one-dimension [4^] . 

Integration of the scaling function with respect to x 
reproduces the scaled velocity distribution &(z). The 
complementary scaled size distribution cannot be found 
in a closed elementary form, except for the special case of 
asymptotically flat distributions (fi = 0) where both of 
the single variable scaling functions are purely Gaussian. 



C. Generalizations 

A natural generalization is to spatially heterogeneous 
initial velocity distributions, Pq(x,v). The time and 
space dependent cluster velocity distribution, P(x,v,t), 
follows from a straightforward generalization of the basic 
derivation in the homogeneous case 



P(x, v, t) = Pq(x — vt, v) exp 



N(t) 



(10) 



fX~V t 

dv' I dx'P (x',v' 
(9) 



For instance, consider the special case where cars are uni- 
formly distributed in the region x < while the region 
ahead is empty. Here, one finds a governing length scale 
x ~ vt ~ t a with the same exponent a as in Eq. (^|) . This 
length scale characterizes a propagating front of clusters, 
and the space dependent concentration c{x, t) becomes 
a function of the scaling variable X = xt~ a , namely 



This growth law is universal as the dependence on the 
details of the initial velocity distribution is secondary, 
entering only via the prefactor. 

In summary, a scaling asymptotic behavior character- 
izes the kinetics of clustering in no-passing zones of one 
lane roadways. The corresponding scaling exponents and 
scaling functions are characterized by the small-velocity 
statistics of the initial velocity distributions. Remark- 
ably, it is possible to derive the exact Boltzmann equa- 
tion in this case. 



III. PASSING ZONES 

We now describe the complementary case of passing 
zones where fast cars can pass slow cars. The model 
we consider is a straightforward generalization of the no- 
passing case. The initial conditions are identical: cars are 
distributed randomly in space with concentration Cq and 
their velocity is drawn from the intrinsic velocity distri- 
bution Pq (v) . The characteristic velocity scale is taken to 
be Vo- In the absence of other cars, cars drive ballistically 
with their intrinsic velocity. In the presence of other cars, 
two competing mechanism may cause a change in the car 
velocity. Collisions lead to slowing down: when a cluster 
overtakes a slower cluster, a larger cluster moving with 
the smaller of the two velocities forms. Passing leads to 
a velocity increase: every car inside a cluster may spon- 
taneously pass the lead car and resume driving with its 
intrinsic velocity. The corresponding passing rate equals 
a constant, . This is a significant simplification; in 
realistic situations only the first few trailing cars may be 
able to pass. 

It proves convenient to introduce dimensionless veloc- 
ity, space, and time variables: v/vq — > u, xcq — > x, 



c v t - 
where 



t. This rescales the passing rate, t 



R 



t, 



c v t 



col 



R- 



(11) 



is the ratio of the passing time t pas = t to the collision 
time t co i = (cqVo)^ 1 . We term this fundamental dimen- 
sionless quantity the "collision number" and denote it R 
as it is reminiscent of the Reynolds number - the small 
R limit is straightforward and the large R limit is char- 
acterized by boundary layers. 

The starting point for kinetic theory is again the clus- 
ter velocity distribution P(v, t). The approximate Boltz- 
mann equation reads 
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dP(v,t) 
di 



= R- 1 [P (v)-P(v,t)} 

- P{v,t) / dv'(v-v')P{v',t). 
Jo 



(12) 



This evolution equation assumes molecular chaos, namely 
that the stochastic passing events effectively mix the ve- 
locities, and therefore, spatial correlations can be ne- 
glected. This is clearly an approximation, as the colli- 
sion integral does not coincide with the exact collision 
integral (|2|) derived in the no-passing case (the R — > oo 
limit). Still, this term reflects the fact that collisions oc- 
cur only with slower clusters and that the collision rate is 
proportional to the velocity difference. The passing term 
is exact since the concentration of slowed down cars with 
intrinsic velocity v equals Pq(v) — P(v,t). 

In contrast with the no-passing case, the process is now 
stochastic in nature and the system approaches a nontriv- 
ial steady state. Setting the time derivative in Eq. ( |l2| ) to 
zero we see that the steady state cluster velocity distri- 
bution P(v) = P(v, t = oo) satisfies the integral equation 



P(v) 



1 + R dv'(v - v')P{v') 



Po(v). (13) 



Given the intrinsic velocity distribution this relation gives 
the final cluster velocity distribution only implicitly. In 
contrast, the inverse problem is simpler as knowledge of 
the final distribution, the observed quantity in real traf- 
fic flows, gives explicitly the intrinsic distribution. We 
confirm that in the limit R — > oo, all clusters move with 
the minimal velocity, while in the limit R — > 0, all cars 
move with their intrinsic velocity P(v) — > Pq(v). 

The integral equation ( p"3| ) can be transformed into 
a differential one using the auxiliary function Q(v) — 
R- 1 + dv'(v - v')P(v'), from which P(v) = Q"{v). 
Thus Eq. ( |l3| ) becomes 



Q(v)Q"(v) = R-'Poiv). 



(14) 



The boundary conditions are Q(Q) = R^ 1 and Q'(Q) = 0. 
The auxiliary function Q(v) gives a comprehensive de- 
scription of the steady state. Calculation of impor- 
tant quantities such as the flux J requires knowledge of 
G(v), the car velocity distribution. This quantity satis- 
fies 1 = JdvG(v) and J = JdvvG(v). Following an in- 
volved calculation that requires solution of a higher order 
velocity distribution p3[ ], the car velocity distribution is 
derived explicitly in terms of Q(v) 



G(v) = P(v) 



1 + R 



dw Pq (w) 



du 



[RQ{u) 



(15) 



Hence, for arbitrary intrinsic velocity distributions, the 
entire steady state problem is reduced to the nonlin- 
ear second order differential equation (p~4|) . Given Q(v), 
steady state distributions such as P(v) and G(v) can be 
calculated using the explicit formulas above. 

Except for a few special cases, one can not solve the 
differential equation (nj) analytically. Nevertheless, the 



formal solution above can be used to evaluate generic 
features of the flow. The dimensionless collision num- 
ber R is extremely useful. For low collision numbers, a 
perturbation solution in powers of R can be constructed, 
as the steady state differs weakly from the initial state. 
For high collision numbers, a boundary layer analysis is 
possible as sufficiently small velocities are not affected by 
collisions. These two limits are quantitatively analyzed 
as follows. 



A. Low Collision Numbers 



The flow characteristics in the collision-controlled 
regime, JScl, can be analyzed systematically as a per- 
turbation series in R. For example, the cluster velocity 
distribution and the car velocity distribution read 



P(v)=P (v) 
G(v) <* P (v) 



1 - R / dv'(v - v')P {v') 
Jo 

poo 

1 + R dv'(v' - v)P (v') 
Jo 



(16) 



Consequently, average quantities such as the flux and the 
average cluster size vary linearly in R in this free flow 
regime, J = Jo — const x R, and (m) = 1 + const x R. 
The proportionality constant in the case of the flux equals 
the variance in the initial velocity distribution, indicating 
that the larger the initial velocity fluctuations, the larger 
the reduction in the flux. 

Therefore, weakly interacting "laminar" flows arise in 
the R — > limit. Technically, the steady state remains 
close to the initial state and a perturbation series in the 
collision number is possible. Here, the assumptions made 
in our theory are justified, as the cluster sizes are small, 
and at the leading order, a simplified model were all cars 
in the cluster can pass coincides with a more realistic 
model where only the first few cars may pass. In fact, a 
basic prediction of the model, namely linear growth of the 
average cluster size with the flux is consistent with em- 
pirical data, obtained from observations of traffic flows in 
a secondary rural road in Los Alamos, New Mexico p3|. 



B. High Collision Numbers 

The limit of high collision numbers corresponds to 
dense, congested flows where large clusters form. The 
analysis in this passing-controlled regime is more sub- 
tle since the condition R L dv'(v — v')Po(v') <gc 1 is sat- 
isfied only for small velocities. No matter how large 
R is, sufficiently slow cars are not affected by colli- 
sions, and P(v) is still given by Eq. ( |l6| ) when v -C v*. 
The threshold velocity v* = v*(R) is estimated from 

R Jo* dv(v* -v)P (v) ~ 1. 
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In the limit R — > oo limit, statistics of the slowest cars 
dominate the flow. Again, it is useful to consider intrinsic 
distributions with an algebraic small velocity form 
For such distributions, the threshold velocity decreases 
with growing R according to 



R f+ 2 . 



(17) 



One can show that the flux is proportional to this 
velocity J ~ v* . For v ^> v* , the collision inte- 
gral in Eq. ( p"3| ) dominates over the constant factor 
and RP(v) f Q dv'(v — v')P(v') ~ . Anticipating an al- 
gebraic behavior for the cluster velocity distribution, 
P(v) ~ R a v s when v 3> v* , gives different answers de- 
pendent on whether the cutoff exponent /i is positive or 
negative. The leading behavior for v 3> v* can be sum- 
marized as follows 



^Yiv/v*)^ 1 /i<0; 
P(v) ~ { (v/v^llniv/v*)}-? fj, = 0; 



(18) 



On the other hand, P(v) = Po(v) for v <C v*. This 
shows that the velocity distribution develops a boundary 
layer structure, the size of which vanishes in the infinite 
collision number limit. Inside the boundary layer, the 
velocity distribution is only marginally lower than its ini- 
tial values, while the bulk of the velocities are strongly 
suppressed. Similar to the threshold velocity v*, macro- 
scopic characteristics of the flow depend algebraically on 
R. For example, the average cluster size is 



R (n+l)/( M +2) fi < . 

{R/lnRf/ 2 fi = 0; 
R 1 ' 2 M > 0. 



(19) 



Two distinct regimes of behavior emerge. For /i > 0, car- 
cluster collisions dominate while for fi < cluster-cluster 
collisions dominate. Interestingly, in the cluster-cluster 
dominated regime, (m) ~ R a with the scaling exponent 
a = (fj, + l)/(fx + 2) as in the no-passing case (||). Thus 
in the passing case the cutoff exponent /i also plays an 
important role in characterizing the behavior. Moreover, 
the steady state behavior is much richer than that found 
for the clustering kinetics. 

Despite the simplifying assumptions, the model results 
in realistic behavior. The overall picture is both familiar 
and intuitive: due to the presence of slower cars, clusters 
form and the overall flux is reduced. For heavy traffic, 
the characteristics of the flow are solely determined by 
the distribution of slow cars. A single dimcnsionless pa- 
rameter, the collision number R, ultimately determines 
the nature of the steady state. 



IV. THE MAXWELL MODEL 

While a comprehensive analysis of the steady state ve- 
locity distributions is possible using the approximate ki- 
netic theory (f||), other important questions such as the 



relaxation toward the steady state and the nature of the 
cluster size distribution [ fl3[ remain unanswered. To ad- 
dress these issues we consider a further approximation 
where the collision rate is taken to be uniform p5| , ^6| . 
This approximation, known as the Maxwell model, is very 
useful in kinetic theory [Q and it has been recently ap- 
plied to granular gases as well p5|-f47j|. In our case, it al- 
lows for a complete exact solution of the time dependent 
behavior, and additionally, it leads to closed evolution 
equations for the cluster-size distribution. 



A. Relaxation 

In the Maxwell approximation, the collision rate v — v' 
in the Boltzmann equation ( ]l2| ) is replaced by a constant 
factor which we set equal to unity. The corresponding 
rate equation for the cluster velocity distribution reads 



dP(v,t) 
dt 



= R- 1 [P (v)-P(v,t)} 

- P{v,t) / dv'P{v',t). 
Jo 



(20) 



Again, the analysis is performed via a properly defined 
auxiliary function, Q(v,t) = J Q dv' P(v' ,t). The con- 
stant collision rate results in simpler differential equa- 
tions, that are only first order in the velocity. The ana- 
log of Eq. (14) is the integrable steady state equation 
Q(v)Q'(v) = R~ 1 Pq(v). The resulting steady state prop- 
erties are governed by R, with a boundary layer structure 
in the large R regime. The quantitative characteristics 
are somewhat different and for example the threshold ve- 
locity decays with R according to v* ~ rather 
than Eq. (jl7|). If, however, the collision rate is properly 
chosen, namely set equal to (v) rather than unity, we 
recover Eq. (p7|). 

Furthermore, the complete time dependence can be ob- 
tained analytically by integrating the partial differential 
equation Q t — R~ 1 Q V — QQ V . In general, the relaxation 
is exponential P(v, t) — P(v, t = oo) ~ /(w)e~*/ r ^- ) , with 
t{v) = R[l + 2RI a (v)r 1/2 where I Q (v) = £dv'P Q (v'). 
The relaxation time depends on the velocity and the col- 
lision number according to 



J R v < v* 

{[R/Ioiv)} 1 / 2 v^v* 



(21) 



While small velocities are governed by practically fixed 
relaxation times, large velocities are characterized by ve- 
locity dependent decay rates. Furthermore, a large range 
of relaxation scales exists, R 1 / 2 < r < R, with larger 
scales corresponding to smaller velocities. Further anal- 
ysis shows that the same relaxation times underlie the 
car velocity distribution. We expect that while the pre- 
dictions of the Maxwell model are only approximate, it 
correctly predicts the existence of a spectrum of relax- 
ation time scales, and that the qualitative nature of the 
time dependent behavior generally holds. 
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B. The size distribution 



V. DISCUSSION 



The size distribution obeys closed evolution equations 
in the Maxwell model and can be solved exactly p5jp6| . 
It can also be used to address the nature of the passing 
mechanism. To demonstrate this we consider the model 
where only the next-to-leading car in the cluster may pass 
and resume driving with its intrinsic velocity. From nu- 
merical simulations of this model, we find two distinct 
phases. In the laminar regime, clusters are generally 
small, specifically the cluster size distribution is expo- 
nentially suppressed for sufficiently large sizes. When the 
collision number exceeds a certain threshold, an infinite 
cluster is formed, i.e., a finite fraction of the cars in the 
system are in the cluster behind the slowest car. Further- 
more, in this jammed phase the size distribution of finite 
clusters has a fat tail close to a power-law, P m ~ m _r , 
with t « 2. 

In the Maxwell model framework, the cluster size dis- 
tribution P m (t) obeys a closed system of rate equations 



dP„ 
~dt 



= R- 1 [P m+1 -P m ]-cP m + ]- PiPj, (22) 



i+j=m 



HP 

-±=R- 1 [P 2 -P 1+c ]-cP 1 . 
at 



(23) 



These equations were derived by enumerating all possi- 
ble ways in which clusters evolve. For instance, consider 
Eq. (|2^). Collisions reduce the density of single cars, and 
the collision rate is clearly equal to c, as it is velocity- 
independent. The escape term in Eq. (|2^) is understood 
by observing that the rate of return of single cars into 
the system is 2P 2 + Y,j>3 p j = P% - P\ + c. Here P 2 
is singled out since passing transforms it into two single 
cars while an escape from larger clusters produces only 
one freely moving car. 

Similar equations were previously studied in the con- 



text of aggregation-fragmentation processes |48 
lizing the approach of Ref. [^8| we find that a p 
sition occurs at R c = 1 
size distribution is 



. Uti- 
rase tran- 

26[. For large m, the steady state 



Pn 



- 3 / 2 [1-(1-R) 2 ] m R<1, 



-5/2 



R> 1. 



(24) 



Hence in the laminar regime, the size distribution decays 
exponentially in the large size limit. In the congested 
phase, the size distribution has a power law tail, and in 
addition there is an infinite cluster that contains the fol- 
lowing finite fraction of cars in the system: 



/ = 



0, R<1; 
1-Er 1 , R>l. 



(25) 



Interestingly, this phase transition is similar to phase 
transitions in driven diffusive systems without passing 
p8| p2| . Furthermore, The formation of an infinite clus- 
ter is reminiscent of Bose-Einstein condensation M,f49l . 



The most important question raised by the above re- 
sults concerns the validity of the "mean-field" Boltzmann 
equation jl^). Although passing is a stochastic mixing 
mechanism that diminishes correlations between the ve- 
locities and the positions of the cars, such correlations 
do exist, and it will be interesting to determine whether 
quantitative predictions such as the scaling behaviors 
( |l7f ) and (|l^) are altered by spatial correlations. Sim- 
ilarly, the collision term in Eq.(2^) is written in a mean- 
field spirit and that may be the reason for the discrepancy 
between the theoretical prediction r = 5/2 and the nu- 
merically observed value r « 2 of the decay exponent 
P m ~ mT T . 

The primary feature of our model is quenched disor- 
der, which manifests itself in the random assignment of 
intrinsic velocities. Road conditions (construction zones, 
intermittent passing zones, turns, hills, etc.) present an- 
other source of quenched randomness in actual roads [^l) , 
which is ignored in our model. Quenched disorder signif- 
icantly affects characteristics of many-particle systems, 
especially in low spatial dimensions p2| . We have seen 
that this general conclusion clearly applies to our one- 
dimensional traffic model. Little is known analytically 
on the influence of the spatial disorder. 

Finally, one may modify the passing rule so that when 
a car overtakes a slow car, it acquires a new velocity 
drawn from the distribution Pq(v) rather than a pre- 
assigned velocity [p3| . This elementary zero-memory 
model remains highly non-trivial even in the collision- 
controlled limit R — > where clustering can be disre- 
garded. The fate of the system is again determined by 
the behavior of the intrinsic velocity distribution near its 
lower cutoff. If Po(v) vanishes in this limit, the system 
reaches a steady state, otherwise, the system evolves in- 
definitely. Specifically, for intrinsic distributions with an 
algebraic small velocity tail (|]) the long-time asymptotics 
of the average velocity reads 



M*)> 




V > 0; 
M = 0; 
-1 < fi < 0. 



(26) 



These results were derived in a simplified Boltzmann 
framework. In particular, the most interesting behav- 
ior in the evolving regime was obtained by assuming 
that as t — > oo, cars can be divided into two groups, 
the small group of "active" cars which move with veloc- 
ities v ~ I and the vast majority of "creeping" cars that 
hardly move at all. We then ignored collisions between 
creeping cars (since their relative velocity is very small) 
and collisions between active cars (since their density is 
small). Thence, the velocity distribution of active cars 
obeys a linear Boltzmann-Lorentz equation which was 
solved to give (|26|). Comparison with results of molec- 
ular dynamics simulations suggests that the mean-field 



G 



theory description is asymptotically exact. It will be in- 
teresting to confirm this result rigorously. 
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